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Abstract 

Context. An important goal within the quest for detecting an Earth-like extrasolar planet, will be to identify atmospheric gaseous 

bio- signatures. 

Aims. Observations of the light transmitted through the Earth's atmosphere, as for an extrasolar planet, will be the first important step 

for future comparisons. We have completed observations of the Earth during a lunar eclipse, a unique situation similar to that of a 

transiting planet. We aim at showing what species could be detected in its atmosphere at optical wavelengths, where a lot of photons 

are available in the masked stellar light. 

Methods. We present observations of the 2008 August 16 Moon eclipse performed with the SOPHIE spectrograph at the Observatoire 

de Haute-Provence (France). Locating the spectrograph's fibers in the penumbra of the eclipse, the Moon irradiance is then a mix of 

direct, unabsorbed Sun light and solar light that has passed through the Earth's atmosphere. This mixture essentially reproduces what 

is recorded during the transit of an extrasolar planet. 

Results. We report here the clear detection of several Earth atmospheric compounds in the transmission spectra, such as ozone, 

molecular oxygen, and neutral sodium as well as molecular nitrogen and oxygen through the Rayleigh signature. Moreover, we 

present a method that allows us to derive the thickness of the atmosphere versus the wavelength for penumbra eclipse observations. 

We quantitatively evaluate the altitude at which the atmosphere becomes transparent for important species like molecular oxygen and 

ozone, two species thought to be tightly linked to the presence of life. 

Conclusions. The molecular detections presented here are an encouraging first attempt, necessary to better prepare for the future 

of extremely-large telescopes and transiting Earth-like planets. Instruments like SOPHIE will be mandatory when characterizing the 

atmospheres of transiting Earth-like planets from the ground and searching for bio-marker signatures. 
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1. Introduction 

Soon the CoRoT (a space program operated by the French Space 
Agency, CNES) and Kepler (a NASA spacecraft searching for 
habitable planets) missions may discover a transiting Earth-like 
planet. It will then be of prime importance to study such an ex- 
oplanet's atmospheric composition, in order to define its evo- 



* Detailed observations as shown in Figs. 9, 10, 11 and 

12 are only available in electronic form at the CDS via 

anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http : //cdsweb . u-strasbg . fr/cgi-bin/qcat? J/A+A/ 



lutionary status as compared to the known Earth one. Through 
these observations the impact of galactic environmental pertur- 
bations as, e.g. the passing through a dense interstellar cloud 
(Vidal-Madjar et al. 1978| l, may be directly evaluated. However, 
the ultimate goal of these studies will be to know if life could 
have emerged elsewhere as it did once on Earth. 

It has been demonstrated, as for the exoplanet HD209458b, 
that it is possible to detect the presence of some species 
in the planetary atmosphere, like sodium, hydrogen, oxy- 
gen, and carbon JCharbonneau et aL2 002 Vidal-Madjar et al. 
2003, 2004). More recently, additional detections have been re- 



Vidal-Madjar et al.: The Earth as a transiting planet 



ported, including HI from recombination via the Balmer jump 
JBallester et al. 2QQT\ . H2 the main gaseous content via Rayleigh 
scattering ( Sing et al. 2008a , 2008b, Lecavelier des Etangs et al. 
2008b), and TiO/VO (Desert et al. 2008). Th e possible signa- 
ture from H2O in HD209458b (IBarman 20071 1 is now strongly 
questioned by new HST (the Hubble Space Telescope) and 
Spitzer (a Space Telescope, studying the universe in infrared) 
observations, revealing on the contrary the presence of high- 
altitude haze as well as CO in the atmosphere of HD 189733b 
jEhrenreich et al. 20071 Lecavelier des Etangs et al. 2008a, Sing 
et al. 2009, Desert et al. 2009). New transiting planets are dis- 
covered from space by CoRot, launched at the end of 2006, and 
by Kepler launched in March 2009. More specifically, space ob- 
servations are reaching high enough accuracy to enable the de- 
tection of Earth-size transiting planets JBorde et al. 20031 Rouan 
et al. 2009, Leger et al. 2009, Queloz et al. 2009) , some of which 
are possibly "ocean-planets" (Leger et al. 2004) or even similar 



to the telluric planets of our Solar System, at distances compat- 
ible with a so-called habitable zone. Obviously, the characteri- 
zation of the corresponding atmospheres would be an exciting 
achievement. 

Knowing that the probability to have a transiting Earth at 
about one AU from its star is of the order of 0.5%, and that 
the number of available planets rises as the cube of their dis- 
tance from the Sun, Ehrenreich et al. (2006) have shown that 
the number of available transiting planets within 60 pc is simi- 
lar to the number of possible targets available if one looks at the 
planetary emissions within a sphere of about 10 pc around the 
Sun (the current expectations of projects like TPF-I (the NASA, 
Terrestrial Planet Finder via Interferometry) or Darwin (the ES A 
equivalent program). However, the present difficulty of the tran- 
siting approach is the lack of photons (i.e. high S/N) necessary to 
analyze spectra of Earth-like planets when using 2-m class tele- 
scopes. Lecavelier des Etangs & Ehrenreich (2005) have shown 
that many hot Jupiters, hot Neptunes, ocean planets, and other 
low-density planets will be observable with larger ground based 
telescopes, like the future European Extremely-Large Telescope 
(E-ELT). Transmission spectroscopy has akeady proved its po- 
tential to probe atmospheres of extrasolar planets. It could very 
well be the first method which will give access to the atmosphere 
of smaller and cooler transiting planets, hopefully Earth-like. 

Observing in the visible and near UV offers a great ad- 
vantage, because of the strong opacity of some specific nar- 
row line species in that spectral range, and because of the of- 
ten more intense stellar flux at these wavelengths. During a 
transit event, only the upper atmosphere of the planet may 
be seen, thus reducing the many perturbing effects potentially 
present when observing at lower resolution, i.e. at lower alti- 
tudes, like the presence of clouds, continents, oceans etc. as 
shown in lArnold et al. 2002llmolf et al. 2002llArnold 2008l and 
references therein. Indeed, in the core of the spectral lines we 
will sample only the upper atmosphere, where one can expect 
that the mixing of different species has already taken place, thus 
revealing an "average" planetary atmosphere representative of 
overall signatures. 

Following the pioneering work of Sagan et al. (1993) search- 
ing for "life on Earth" from remote sensing, our purpose is 
to observe the Earth as a transiting planet in order to bet- 
ter test the present models of Earth-like planets and thus get 
ready to properly analyze future observations of real transit- 
ing Earth-like extrasolar planets. Using the Moon as a reflec- 
tive surface during total or partial eclipses provides the needed 
situation which furthermore can be exploited from the ground. 
We have conducted a first attempt during the 2008 August 16 



lunar ecUpse with the SOPHIE high-resolution spectrograph 
dPerruchot et al. 20081 [Bouchy et al. 2 009) at the 1.93 cm tele- 
scope of the OHP (Observatorre de Haute Provence) observatory 
covering the 3 900 to 7 000 A spectral range at a resolving power 
of R~75 000. The choice of a high-resolution spectrograph al- 
lows us to better identify absorbing species via their detailed 
spectral signatures. 

Other observations of this eclipse have been conducted 
dPalle et al. 2009l l over a wider wavelength range (0.36- 2.4 /zm), 
but at a lower spectral resolution (R< 1 000), resulting in a 
mainly qualitative identification of the Earth's atmospheric com- 
pounds. Here, we will focus on quantitative results. In Sect. |2] 
the SOPHIE observations are described, then a quantitative data 
analysis is given in Sect. [3] followed by a direct comparison with 
current model predictions in Sect.|4] We show that indeed, in the 
visible range some of the main Earth's atmospheric compounds 
are identified and evaluated. 



2. Observations 

The observations were conducted with the cross-dispersed, 
environmentally stabilized echelle spectrograph SOPHIE, 
dedicated to high-precision radial velocity measurements 
dPerruchot et al. 20081 Bouchy et al. 2009 1. They were secured 
in high-resolution mode, i.e. the spectrograph was fed by a 
40 jt/m slit and an optical scrambler located at the output of the 
optical fiber, allowing a resolving power A/AA = 75 000. 

Two optical fibers, 3 arcsec wide, separated by 1.86 arcmin 
were used. Both were placed on the Moon and aligned along the 
east-west direction. 

The spectra were extracted from the detector images with the 
SOPHIE pipeline, which includes localization of the orders on 
the 2D-images, optimal order extraction, cosmic-ray rejection, 
wavelength calibration, corrections of flat-field and charge trans- 
fer inefficiency at low signal-to-noise ratio (Bouchy et al. 2009)l. 



2. 1 . Fibers position on fhe Moon 

Table [T] lists the three main spectra used in our analysis along 
with their positions over the Lunar disk as shown in Fig. [1] Only 
one observation during the eclipse was useful, taken at 20h44 UT 
in the penumbra, close to the Earth shadow. Figure [1] shows 
where the fibers were positioned, i.e. over high-albedo Imbrium- 
type or pre-Imbrium-type terrains (for more information see the 
URL at the end of the reference list), which maximized the flux 
during the eclipse. The two other observations listed in Table [1] 
are two among a total of 1 3 calibration spectra recorded after the 
eclipse (full Moon outside the penumbra) at similar air masses 
(sec(z) in Table [1] where z is the zenith angle) and over the two 
different regions pointed to during the eclipse. Unfortunately, af- 
ter the eclipse the guiding camera was saturated and an accurate 
position of the fibers could not be checked and adjusted visually 
as during the eclipse. Nevertheless the telescope coordinates in- 
dicate that the fibers were positioned over the same kind of lunar 
terrains. 



2.2. Observations log 

We collected only sparse observations, because clouds appeared 
just after the first exposure taken at 20h44 UT. These clouds dis- 
appeared for a few minutes at 22h44 UT, just before the umbra 
left the Moon, though these observations are more noisy and are 
still possibly affected by intervening clouds. 
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Table 1. Fibers positions on the Moon during and after the eclipse. Times are for mid-exposure. 

Date; Fiber A location Fiber B location Exposure BERV Air mass 

Time (UT) longitude; latitude longitude; latitude time (s) (km/s) sec{z) 

2008-08-16; 20:44:21 
2008-08-17; 00:11:56 
2008-08-17; 02:55:45 



E 52.0°; N 9.8° 


E68.5°;N 12.1° 


198.92 


-0.064 


2.94 


E 39°; 8 51° 


E50°;S48° 


2.00 


-0.408 


1.86 


E47°;N20° 


E 64°; N 22° 


3.00 


-0.667 


2.99 




Figure 1. Fibers A and B on the Moon in penumbra. The umbra 
is the dark area at the left (upper figure at 20h44 UT) and at right 
(lower figure at 22h44 UT). The distance between the fibers is 
1.86 arcmin, east- west oriented. 



The observations completed at 20h44 UT were done with 
fibers A and B at about 2 and 4 arcmin from the penumbra/umbra 
limit respectively. 

Spectra were also taken later in the night after the eclipse, 
when the clouds moved away. In particular at 02h56 UT spectra 
were taken over the same location on the Moon and at the same 
air mass as those at 20h44 UT. The spectrum taken at 00hl2 UT, 
although pointed over another region of the Moon, is also listed 
because it was taken with a shorter time gap from the 20h44 UT 
eclipse observation, leading to a smaller differential BERV cor- 
rection ("Barycentric Earth Radial Velocity" correction of the 
Earth's motion along its orbit to evaluate the spectral observa- 
tions in heliocentric coordinates as automatically done through 
the SOPHIE pipeline). The corresponding BERV shifts are also 
listed in Table [T] 

The calibration spectra (taken later in the night after the 
eclipse) are used to correct the eclipse spectra from the direct at- 



mospheric absorption along the line of sight between the Moon 
and the OHP Observatory and to extract the altitude at which the 
Earth atmosphere becomes transparent enough for a grazing line 
of sight. 

3. Data analysis 

Because we are searching for atmospheric signatures in the gath- 
ered spectra, we have first to recalculate the wavelengths in the 
geocentric reference frame. We know that there all atmospheric 
signatures will be present at identical wavelengths, while the so- 
lar line signatures will be slightly displaced by an amount re- 
lated to the spectral shift induced by the BERV. Note however 
that at the SOPHIE resolution of about ~4km/s, these shifts of 
less than 0.5km/s are likely insignificant. The possible impact of 
the BERV effect on our data is assessed in Sect. |5] 

Our observations present two types of Earth atmospheric ab- 
sorption signatures : i) the "vertical" ones produced along the 
path length between the Moon and the telescope and ii) the "hor- 
izontal" ones due to the solar light grazing the Earth atmosphere 
before reaching the Moon. We are only interested in the latter 
and will therefore have to correct the data from the "vertical" 
signatures. 

Finally, knowing that the SOPHIE pipeline is optimized for 
high radial velocity precisions and not for flux evaluations, we 
can expect some difficulties related to flux estimates, i.e. echelle 
order corrections as well as precise zero level evaluations. 

3.1. The "vertical" Earth atmospheric transmission function 

The "vertical" atmospheric transmission is a function that de- 
pends upon the air mass AM, noted Tam(^)- Along the line of 
sight, an initial spectrum I{A) reaching the top of the atmosphere 
is transformed into an observed spectrum OamW ■ 



OamW = 1(A) X TamW- 



(1) 



For a plane parallel atmosphere (an approximation valid to air 
masses smaller than ~5) the air mass AM is equal to a=sec(z) 
and the atmospheric transmission function is then simply T''{A), 
the transmission function T at the power a ( IBird et al. 1982l l. 
where T is the transmission in the vertical direction (air mass 
of 1). 

Indeed for an air mass a equal to the sum of two air masses ai 
and fl2, one can evaluate the observed spectrum in the following 
manner : 



C>fl,+„, = / X T'"^"'- = (/ X r"') X T"'-, 



(2) 



showing that the transmission function T(A) follows the condi- 
tions : 



r'"+"2(i) = T'"(A)xT"'(A) 

and 

T^\A) = 1. 



(3) 



(4) 
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If the transmission function T(A) is known at a given air 
mass, it can be evaluated at any other air mass. Noting that 



0„,(A) _ T"'(A) 



_ Y'Ci 



"HA), 



(5) 



we were able to evaluate directly the transmission function at 
different air masses through the complete series of spectra gath- 
ered during the full Moon observations. We normalized all eval- 
uations extracted from observational ratios to an air mass equal 
to 1 by calculating 






l/(ai-fl2) 



= (r«l-«2(^)) 



1/(01-02) 



= T(A). 



(6) 



From the 13 full Moon observations OJA), we extracted 24 
couples well separated in air mass each containing independent 
evaluations of the transmission function T(A). In Fig.|2] the av- 
erage of all evaluations are shown, providing the reference trans- 
mission function T(A) that we will now use in the following 
data analysis to re-evaluate all observed eclipse spectra as if they 
were observed from outside the atmosphere. 




WAVELENGTH (A) 

Figure 2. Average atmospheric transmission function T{A) as 
evaluated during full Moon just after the eclipse supposed to be 
similar to the one during the eclipse observations. The dashed 
line shows a transmission model ( |Hayes and Latham 1975 1, 
which includes both Rayleigh diffusion and aerosols (optical 
depth of 0.035 at 5 320 A typical for clear OHP nights). The 
experimental evaluation was normalized in order to match the 
model at our reference wavelength Aq - 4530 A (see Sect. [S]). 
O2 as well as H2O absorption bands are clearly seen. Note that 
the match with the model is good except for the central region 
from 5 000 to 6 700 A, where a clear additional absorption is de- 
tected. This is the Chappuis band from ozone (noted O3), not 
included in the model calculation. 



3.2. Other corrections of the SOPHIE spectra 

One can note in Fig. |2] the SOPHIE echelle order signatures as 
regular wiggles over ~80A from 4000 to 4500A, where they are 
particularly visible. Those instrumental signatures should not be 
seen in an atmospheric transmission function. These will cer- 
tainly be a perturbation in our data analysis, but for the time 




5892 

WAVELENGTH (A) 



Figure 3. At 20h44 UT, the (Fiber B / Fiber A) ratio (solid line) 
is shown in the Nal doublet spectral region. The Fiber A spec- 
trum (dashed line) with its corresponding zero level (horizontal 
dashed line) as extracted from the SOPHIE pipeline as well as 
an H2O model (dotted line) are overplotted (in arbitrary scales 
simply to fit within the figure frame) to show the positions of the 
solar as well as H2O spectral lines in the spectrum. Most of the 
observed lines are due to atmospheric water vapor. They disap- 
pear in the fiber's spectral ratio, while the broad spectral wings 
of the two Nal solar lines are still clearly visible. 



being we will try to analyze the data without a modification to 
the SOPHIE pipeline. We will however evaluate the impact of 
this approximation a posteriori by looking how these order sig- 
natures are affecting the results (see Sect.|5]l. 

In Fig. [3] both the direct Fiber A spectrum observed dur- 
ing the eclipse as well as the Fiber B/Fiber A flux ratio are 
shown close to the Nal doublet solar lines. Note that the numer- 
ous spectral lines, mostly due to H2O in the Earth atmosphere, 
completely disappear from the Fiber B/Fiber A ratio. This is not 
surprising because both Fibers A and B are observed through the 
same atmospheric (and thus air mass) layer. 

Although most of the weaker solar lines are also erased in 
this ratio (as it can be directly seen in Fig. |2l which also shows 
spectral ratios, see Eq. |6ll, for strong and deep solar lines like 
e.g. the Nal lines, as shown in Fig. [3] the signature of their broad 
spectral wings is still visible in the Fiber B/Fiber A ratio. This 
should not be the case if Fiber A and B have simply a different 
transmission coefficients, because they are both simultaneously 
observing (almost) the same solar spectrum. 

The Nal lines (reaching close to the zero level, see Fig. [3]) 
are more sensitive to any unaccounted for, even small, zero level 
shift. This is not the case for the weaker solar lines. 

Knowing that a zero level shift is present in the extracted 
spectra, we search below for some plausible explanations. 

3.2.1. An instrumental zero level shift 

Over the SOPHIE CCDs, the lowest counts are not equal to zero 
(even where no counts should be registered), as for instance at 
the shortest wavelengths, where the instrument sensitivity drops 
sharply. For this reason, we assumed the possibility of a zero 
offset in the SOPHIE extracted counts. We thus subtracted a 
value from the SOPHIE spectra to make the broad Nal solar lines 
wings disappear (as they should). This correction is compatible 
with the data, as reported in Table |2l 
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Table 2. SOPHIE fluxes as extracted from the data sets, both at 
the minimum of deep absorption signatures and where the flux is 
maximum. Times are for mid-exposures. Flux units are arbitrary 
but comparable in relative terms because they correspond to the 
raw data (CCD ADU) corrected for by the same pipe line. 



Fiber ; 
time (UT) 


Min. flux 
inNal 


Min. flux 
in02 


Max. 
flux 


A ; 20:44 


0.0007 


0.0005 


0.0205 


B ; 20:44 


0.004 


0.002 


0.106 


A ; 02:56 


0.0018 


0.0010 


0.0514 


B ; 02:56 


0.0016 


0.0008 


0.0471 



We tried three zero shift level scenarios : 

1 . we assumed that stray light due to the nearby full Moon is 
present. However, subtracting various proportions of the full 
moon spectrum to the other spectra did not induce any atten- 
uation of the Nal solar wings and on the contrary produced 
many additional solar lines; 

2. we supposed that a constant value has to be subtracted from 
all data, proportional to the total flux gathered in each spec- 
trum. Again no satisfactory solution was found; 

3. we assumed that a unique constant value had to be used for 
all spectra, a correction only related to the detector and in- 
dependent of the total flux. Then we found that by subtract- 
ing 0.0003+0.0001 from afl SOPHIE extracted values, it was 
possible to completely wash out the Nal solar wings. This 
0.0003 value represents a maximum relative continuum cor- 
rection of 14% in the blue, decreasing to 1.5% in the red, for 
the fiber A spectrum. For fiber B, which is brighter, the rela- 
tive continuum correction is one order of magnitude lower 

We verified that this last correction (3.) was also adapted at 
the end of our study, because no broad wing solar Nal lines sig- 
natures should show up in any information related uniquely to 
the Earth's atmosphere. The corresponding final result (with and 
without this correction) is shown in Sect.|5] Fig.fTOl 

This simple third correction (3.) does not seem to be 
wavelength-dependent, because as we show in Table |2l at the 
wavelengths of the very deep O2 absorption lines near 6 900 A, 
it is clear that this correction could not be larger than the min- 
imum flux seen at the bottom of the absorption lines, and thus 
has to be somewhere in the to 0.0005 range, which we found 
to be the case. A posteriori, we will also see that the final results 
are all compatible with this correction. 

After empirically evaluating the type of correction needed, 
we found that another possible explanation of atmospheric origin 
may be given and we develop it in the next section. 

3.2.2. An astrophysical zero level shift : the ring effect 

Indeed, following a suggestion of our referee, we recalled that 
some of the sunlight going through the Earth atmosphere is re- 
fracted in the lower layers of the atmosphere, producing for an 
observer on the Moon within the Earth umbra, an emission ring 
all along the Earth limb. 

This ring effect can contribute here in two ways. First, it can 
add a constant intensity (a few percent shift of the zero level) 
to sunlight as seen in scattered light from the sky, to the ex- 
tent that the penumbral light contains not just transmitted light, 
but also forward-scattered light ; this should slightly dilute the 
penumbral spectrum and introduce a zero level correction of type 



(2). Second, the Ught reflected from the lunar surface itself will 
contain the ring effect, again a few percent additive background 
spectrum, a correction of type (2.) or (3.). 

The ring relative intensity and spectral signature were al- 
ready observed during Lunar eclipses through the study of the 
lunar light reflected from the umbral regions, only lit by the 
Earth ring. Danjon (1936) very early understood that the umbral 
luminosity was related to the transparency of the Earth lower at- 
mosphere and after observing several lunar eclipses, he proposed 
a related five level scale for the lunar umbral luminosity. 

Observation revealed (as in e.g. Hernitschek et al. 2008) that 
a 5 magnitudes drop exists at the penumbra-umbra (PU) tran- 
sition, showing that the ring contribution to the penumbral flux 
is on the order of 1 %. Furthermore, the ring contribution de- 
creases when moving away from the PU transition deeper within 
the umbra, as the expected signature of forward-scattered light. 
The ring contribution within the penumbral region is indeed pro- 
portional to the length of the Earth's limb (Lg) just in front of the 
solar disk appearing on top of the Earth limb. 

Thus, when moving in the penumbral domain away from 
the PU transition, the direct solar contribution increases like 
the surface of the solar disk above the Earth limb (~L|, see 
Fig.HJl, while the ring contribution increases as the length of the 
Earth's limb covered by the solar disk (~Le). This means that 
for our observations, the relative ring effect contribution should 
be stronger in Fiber A than in Fiber B according to the geome- 
try of our 20h44 UT observation, for which Fiber A, presenting 
a lower signal, is closer to the PU transition than Fiber B. The 
ring effect contribution in the reflected Lunar light seems thus to 
better correspond to our empirical correction scenario (3.). 

This is confirmed in the Gedzelman & Vollmer (2008) sim- 
ulation of irradiance of lunar eclipses. They furthermore show 
that the ring contribution at the PU transition should present an 
almost flat spectral signature from 4000 to 7000A (see their Fig. 
1), i.e. over the whole SOPHIE spectral range. This was directly 
observed during the 2008 August 16 lunar eclipse by Palle et 
al (2009, see top of their Fig. S3 of the supplementary section), 
who detected only weak spectral signatures deeper in the umbra 
as also predicted by Gedzelman & Vollmer (2008). 

In summary the ring contribution should be on the order of 
1 % of the signal in Fiber A, and be at a relatively lower level 
in Fiber B with no wavelength dependence. This is almost ex- 
actly what we have found from our empirical study (3.), which 
presents the same 0.0003 zero level shift in both Fibers, i.e. in- 
deed on the order of 1 % in Fiber A and relatively less in Fiber B 
(see Table |2]i. Additional discussions in Sect. |5] will show that 
this correction is on the same order near the Nal doublet lines at 
5 900 A and near the molecular O2 lines at 6 900 A and thus is 
indeed weakly variable with wavelength. 

The ring effect could also explain why the ratios (Eclipse 
A)/(Full Moon A), (or B) could also be affected, ratios that we 
are going to use and describe in the following sections. 

Both these instrumental and ring effects are probably present 
in our study and need to be corrected for empirically. In future 
studies we will also try to observe umbral spectra in order to 
have a more quantitative evaluation of the ring effect and thus 
include its contribution in our analysis. 



4. Quantitative analysis 

We present in this section the approach developed to extract a 
quantitative information from the observed spectra related to the 
altitude at which the Earth atmosphere becomes transparent as a 
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function of wavelength, a result very similar to the one obtained 
for a transiting planet. 

4.1. Determination of the unabsorbed full Moon and eclipse 
spectra 

These unabsorbed spectra are simply recomputed by dividing the 
observed spectra Oa{^) by the transmission function T{A) at the 
power of the precise air mass of the corresponding observation : 



1(A) = Oa(A)/r(A). 



(7) 



All these calculated corrected spectra are those an observer 
should have obtained from above the Earth atmosphere. We will 
note these spectra E(A) and F(A), corresponding to the "eclipse" 
and "full Moon" corrected spectra respectively. 

4.2. Earth's atmosphere effective thickness 




Figure 4. Geometry of a crescent of Sun, S , above the Earth's 
reference limb (Os and Rs are the center and radius of the Sun, 
Oe and Re those of the Earth). The scene (not to scale) is ob- 
served from the Moon at the location where are positioned fibers 
A or B within the penumbra. From the observational knowledge 
of Sa or S b as extracted from Eq. [TT](see text), the entire de- 
fined geometry of the observations allows a direct evaluation of 
the corresponding arc lengths. La or Lb- 



From the deep penumbra, where the spectrograph's fibers are 
located during the observation, a crescent of Sun is seen above 
the Earth's limb. The Moon irradiance is thus a mix of direct 
Sun light (i.e. unaffected by the Earth's atmosphere), and solar 
light that passed through the atmosphere, mostly along an arc L 
where the solar disk intersects the Earth's limb (Fig.|4|. This mix 
of direct and absorbed solar light essentially reproduces what is 
recorded during a transit. Let us note 5 © as the surface of the full 
solar disk, S the fraction of solar surface visible above Earth's 
geometrical limb during the eclipse, and h is an equivalent height 
over which the atmosphere can be considered opaque, at each 




Figure 5. L - f(S ) function computed here for the geometry 
of the eclipse: S is the surface of the solar crescent above the 
Earth's reference limb - as observed from the Moon at the loca- 
tion of either fiber A or B within the penumbra. From there, the 
diameters of the Earth and the Sun are 6 768.0 and 1 890.6 arcsec 
respectively. This geometry implies the relation shown between 
the length of the Earth's hmb covered by the Sun L, and the 
surface of the solar crescent, normalized to the total solar disk 
surface ^o as seen from the Moon. Since Sa/Sq or S b/Sq are 
known from observations, the corresponding arcs lengths La or 
Lb are also known, and thus h{A) evaluated. 

wavelength, from the point of view of the geometrical cross sec- 
tion. 

The total column density along a grazing line of sight pass- 
ing through the terminator at different altitudes could be esti- 
mated following the Fortney et al. (2005) formulation: the opti- 
cal depth, T, in a line of sight grazing the Earth limb at an altitude 
h is given by 



t(A, h) K o-{A)n(h) ^l2nREH, 



(8) 



where Re is the Earth radius, H the atmosphere scale height, 
and n{h) - n(h=())exp{-h/H) the volume density at the altitude h 
of the main absorbent with a cross section cr{A). The scale height 
is given by the relation H = kT/fug, where k is the Boltzmann 
constant, T the temperature, fi the mean mass of atmospheric 
molecules times the mass of the proton, and g the gravity at the 
Earth radius, g = MeG/R\ with Me equal to the Earth mass and 
G the gravitational constant. 

For any wavelength, the line of sight becomes opaque at 
an effective altitude h for an effective optical thickness of 
Teff = 0.56 as demonstrated to be appropriate for most atmo- 
spheric conditions by Lecavelier des Etangs et al. (2008a). We 
will then be able to compare our equivalent effective altitude h 
as a function of wavelength to model calculations of effective 
altitudes. 

The effective altitude of the atmosphere at a wavelength A 
is calculated by solving the equation t(A, h) = T(,//. Using the 
quantities defined above, the effective altitude h is given by Eq|9] 



h(A) = //In {ing,=o)(TiA)/T,ff) ^|2nRE|kT^lg) . 



(9) 



In principle, a simple relation links the two E{A) and F(A) 
corrected spectra due to the geometry of the problem 



E(A) = F(A) X 



S -Lx h(A) 



(10) 
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4.3. Evaluation ofh(A), the effective tieiglit of tlie Earth 
atmosphere. 




WAVELENGTH (A) 

Figure 6. Absorbing atmosphere thickness versus wavelength, 
evaluated according to Eq. [14] for the 20h44 UT eclipse obser- 
vation associated to the 02h56 UT full Moon ones. The refer- 
ence altitude equal to km (dashed line) has been chosen for 
A - 4530 A. The profile is dominated by the Chappuis band 
of ozone from 4 800 to 6 900 A, while some signatures due to 
sodium and oxygen show up near 5 890 A and 6 900 A respec- 
tively. In the blue region (4000 to 4 800 A), the Rayleigh scat- 
tering is clearly seen. Across that region, the signal drops and 
produces more noise over the extracted altitudes, also revealing 
there that the SOPHIE orders are not fully corrected. This is also 
true longward from about 6 500 A. 

We can define a reference altitude /z(/lo) = at an arbitrary 
wavelength Aq. Then, S is more precisely the fraction of solar 
surface visible above the Earth's limb corresponding to our ref- 
erence altitude /i(/lo)- It follows that Eq.[TO]reduces to 



E{Ao) = F{Aq) X 



(11) 



The S ISq ratio can be geometrically calculated for a given 
observation time, but may also be directly estimated from the 
E{Ao) and F{Ao) values of the flux measured during and outside 
the eclipse at Aq. We decided to take the reference wavelength 
over a 20 A range centered at /io = 4 530 A, where the absorption 
in the Earth atmosphere (see Fig. |2]i is mainly due to Rayleigh 
scattering. Because the altitude above sea level (corresponding 
to the same reference wavelength) is the standard outcome from 
any atmospheric model calculation (see next section), it will be 
then simple to shift all our evaluated altitudes relative to our ref- 
erence altitude h{AQ) into "real" altitudes, above sea level. Indeed 
all our altitude evaluations h{A) relative to h{A()), could be posi- 
tive or negative depending on more or less absorption at A than 
dXAo. 

The ratio of Eqs.[TO]and[TT]gives for Fiber A 



EaJA) Fa{X) i-A .,.,1 

- X [1 X «(/!)], 

Ea{Aq) Fa{A(,) Sa 

and for Fiber B 
Eb{A) Fb(A) 



Eb(Aq) Fb(A()) 



x[l 



X h(A)], 



(12) 



(13) 
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Figure?. Same as Fig.|6] but here the h{A) variations are binned 
over 200 pixels (~2 A) for clarity. Here the altitude variations 
are, following Eq. [14] extracted from the eclipse observations of 
20h44 UT with either the 02h56 UT (soUd line) or the 00hl2 UT 
(dotted line) full Moon ones. Without the binning the difference 
was not visible. Here we can note that both possible perturba- 
tions related to either the BERV shift as well as the air mass 
correction are indeed negligible (see text). On the contrary some 
broadband variations could be seen, probably linked to albedo 
variations of the different lunar regions where the fibers were 
pointed (see Table [D- In the blue region (4000 to 5 000 A), the 
fringes are due to an imperfect correction of the SOPHIE orders 
and not to classical CCD interferential fringes (see text). 

in which La and Lb are the lengths of the arc of limb at the 
h(Ao) level as seen from the locations of Fiber A and Fiber B 
respectively, as projected on the Moon, and similarly Sa and S b 
are the corresponding surfaces of the Sun above these Earth's 
limbs. 

The difference between the ratios for Fibers A and B, (Eq.[T2l 
-Eq.[T3]i, is 



Ea{A) ^^ Fa(Ao) Eb{A) ^^ Fb(Aq) 
Ea(Ao) ^ Fa(A) Eb(Ao) ^ Fb(A) 



Lb La . 

— - — I X h(A). 



(14) 



Therefore, the measurements (left side of Eq.[T4li are propor- 
tional to h. 

However, to fully solve the problem, one needs to evaluate 
the factor of proportionality in front of h. We know the rela- 
tive angular diameters of both the Earth (above the Earth's ref- 
erence limb) and the Sun, as seen from the Moon at the time of 
the observations, which unambiguously defines (geometrically) 
how far above that selected Earth's limb reference the solar disc 
emerges. From a simple geometrical calculation (see Fig.|4|, one 
can evaluate the corresponding arc lengths La and Lb, which are 
directly derived from the relation L = f(S ) that we have calcu- 
lated for the time of the observations, as shown in Fig. [5] Finally, 
the direct measurements of 5 a and 5 b are extracted from the ob- 
servations through Eq. [TT] 

For example, from our best observations made on August 16 
at 20h44 UT during the eclipse, compared to the full Moon ob- 
servations made on August 17 at 02h56 UT, we find the follow- 
ing corrected flux ratios in the 4 520 - 4 540 A reference spectral 
band : 



Vidal-Madjar et al.: The Earth as a transiting planet 



EAi'k))/FA('lo) - 0.44 % a value similar to Sa/Sq as seen 
from the Moon where Fiber A is located; this leads to La - 
569 arcsec; 

E b('^q) / F b{^o) - 3.12 % a value similar to S b/Sq as seen 
from the Moon where Fiber B is located; this leads to Lb - 
1 063 arcsec. 

In conclusion, from all parameters fully known at a given 
observing time, Eq.[T4]gives access to direct evaluations of h(A), 
either in arcsec seen from the Moon or in km above (or below) 
the Earth's reference limb at h(Ai)). 

The extracted values of h(A) are shown in Fig. |6] obtained 
when using the 20h44 UT eclipse observations compared to the 
02h56 UT full Moon ones at about the same air mass and over 
the same regions on the Moon. The extracted altitudes are rela- 
tive to the reference altitude h{Ao), here equal to 0. As expected, 
the evaluations can be either positive (higher altitudes) or nega- 
tive (lower altitudes). 

To check the validity of several of our assumptions we calcu- 
lated the same h{A) function extracted again from the 20h44 UT 
eclipse observations compared to both full Moon observations, 
the one at 02h56 UT as well as the 00hl2 UT one made at dif- 
ferent air masses, BERV and over distant regions on the Moon 
(see TablelTJ. This allows us to check that our evaluations are not 
sensitive to air masses or BERV variations (see Fig.|7]i, because 
both extracted h(A) functions are very similar in all their details. 
Although nearly identical (non discernable at full spectral reso- 
lution as in Fig.|6]l, the result indicates however that some broad 
band variation is present. It is not caused by an improper air 
mass correction because the produced effect cancels in two spec- 
tral regions: the first one at 4530 A (by definition, since this is 
our reference wavelength) and again around 6 300 A. That the 
observed relative variation of the transmission remains small in 
the blue and more important in the red is indeed in contradiction 
with an air mass correction, which has to be stronger in the blue 
(see Fig. |2]). We thus conclude that the observed relative varia- 
tion is more probably due to lunar albedo changes over different 
lunar regions. 

This broadband effect could however be the cause of addi- 
tional systematic errors in the altitude evaluations. 



4.4. Error estimation 

A simple visual analysis of Figs. |6]and|2]gives some idea of the 
errors estimated for the h(A) evaluations. First in Fig. |6l the er- 
ror due to the photon noise in the observations is ~ ±1 km for 
all data points shown but clearly increases toward shorter wave- 
lengths (below 5 000 A) due to the decreasing sensitivity of the 
instrument toward these wavelengths. Furthermore in this spec- 
tral region the impact of the improper order corrections is more 
and more important as more easily visible in Fig. |7] This ef- 
fect induces an additional systematic error below 5 000 A (of 
~ ±2.5 km) and above 6400 A (of ~ +1 km). Finally, although 
not all possible lunar terrains were sampled, we can guess from 
the different regions observed that an additional systematic error 
due to lunar albedo variations could be on the order of ~ + 1 km. 
This means that for our first attempt to estimate the Earth h{A) 
function we cannot produce absolute altitudes to better than ap- 
proximately ~ ±2.5 km. 



5. Discussion 

5.1. Model calculation 

The information extracted from the presented spectra is h{A), 
thus we are exactly in the situation of an extrasolar planetary 
transit. In effect, a spectrum recorded during a transit is not the 
spectrum of the absorbing species, but only the spectral sig- 
natures (in terms of altitude) of the highest absorbing species 
within the observed atmosphere. For instance, in the central re- 
gion of the spectrum (Fig. |6]l the main feature is the Chappuis 
band of ozone, corresponding to altitudes known to be in the 
30 km range, other detectable species are only those able to ef- 
ficiently absorb above. As we will see, oxygen and sodium are 
also detected with the present observations, but not H2O. 

In order to properly interpret our observations, we have to 
use model calculations to estimate the altitude at which the Earth 
atmosphere becomes transparent, as a function of wavelength 
and for a grazing line of sight exactly as in a transit situation. 

Standard Earth atmospheric model calculations are from e.g. 
Ehrenreich et al. (2006) or Kaltenegger & Traub (2009). From 
these model calculations, one can evaluate the effective Earth ra- 
dius at each wavelength, and from it, the effective height of the 
Earth atmosphere as a function of wavelength. From Fig. 3 of 
Kaltenegger & Traub (2009), the effective altitude correspond- 
ing to the tip of the Chappuis band absorption (near 6 000 A), is 
on the order of 30 km, revealing that the atmosphere becomes 
nearly opaque below this altitude at this wavelength. 

In our model calculation we used an unidimensional sin- 
gle scattering transmission model based on the model de- 
scribed by Ehrenreich et al. (2006) to obtain the theoreti- 
cal transmission spectrum of the terrestrial atmosphere. The 
pressure, temperature, and mixing ratio profiles of the consid- 
ered atmospheric components are shown in Fig. [8] and come 
from the US Standard Atmosphere 1976 spring-fall pressure- 
temperature profile dCOSEA 19761 ICox 20 00). For simplicity, 
we chose to only include molecular nitrogen (N2), molecular 
oxygen (O2), and ozone (O3) in the model. Molecular nitro- 
gen contributes to the transmission spectrum via the Rayleigh- 
scattering of light. The Rayleigh-scattering cross section of N2 
depends on the N2 refractive index, which is calculated accord- 
ing to Sneep & Ubachs (2005). Within the SOPHIE spectral 
range, molecular oxygen contributes to the transmission spec- 
trum through i) photoabsorptions by the forbidden ' l,g-^l,„ tran- 
sition bands B (1-0) at 6 880 A and y (2-0) at 6 280 A, and 
ii) through Rayleigh scattering. The line parameters for the O2 
bands are extracted from the HITRAN 2008 molecular spectro- 
scopic database (iRothman et al. 2009l l. scaled to the atmospheric 
temperature and pressure profiles used, and convolved with a 
Gaussian profile with full-width at half-maximum matching the 
spectrograph resolution. The Rayleigh-scattering cross sections 
of O2 are calculated following Sneep & Ubachs (2005) using re- 
fractive indices from Bates (1984). We retrieved the UV/visible 
photoabsorption cross section of O3 at 293 K and 1 bar from the 
GEISA 1997 database (Jac quinet-Hussonet al. I999| l. 

We recall here that we are not doing a fit of our data but just 
compute altitude signatures related to a "standard" Earth atmo- 
spheric model. 

5.2. Broadband signatures 

The variable effective altitude h (as a function of wavelength) 
is shown on Fig.|9]for both model calculations and as extracted 
from the observations, the overall contribution of O3, O2 and 
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Figures. "Standard" atmosphere used in the to model calculation (see text) : Upper left the pressure profile, upper right the 
temperature profile, lower left the N2 (blue), O2 (green) and O3 (pink) mixing ratio as a function of altitude (km) and lower right 
the N2 (blue) Rayleigh, O2 (green) Rayleigh and molecular bands and O3 (pink) cross sections as a function of wavelength (in nm). 



Rayleigh scattering due to both N2 and O2, along with the 
Rayleigh contribution alone are also shown. 

In the selected reference domain from 4 520 to 4 540 A, the 
model calculation predicts an altitude equal to 23.8 km. Because 
our observational evaluation of h is relative to the altitude in 
the reference domain, we shift our relative h{A) evaluations by 
+23.8 km to be comparable to the model calculation. Since fur- 
thermore our reference domain falls at the tip of a SOPHIE 
echelle order fluctuation, we corrected by an additional system- 
atic shift of about 1.5 km (half of the echelle fluctuation). The 
model and observations are shown in Fig.|9] where a 25.3 km 
reference altitude is used in the 4 520 to 4 540 A reference do- 
main range. 

Our observed evaluations as shown in Fig.|9]are binned over 
200 pixels (~2 A) to better underline the observed broadband 
variations. It can be seen that the errors over the broadband sig- 
natures are mainly due to systematic errors, as opposed to pho- 
ton noise, and are on the order of the SOPHIE approximate order 
correction, i.e. ~ ±2.5 km. 

The comparison between the observed variations and model 
calculations is very satisfactory. Even the observed ozone distri- 
bution seems probably less extended in altitude than the standard 
model by about 2 km, a very plausible situation because ozone 
variations are known to exist around the Earth in terms of both 



location and epoch (Millier et al. 1979, Thomason & Taha 2003, 
Borchi & Pommereau 2007). There is no reason that the selected 
standard atmospheric model exactly corresponds to the Earth at- 
mospheric conditions at the precise location over the observed 
Earth's limb. 



5.3. Narrowband spectral signatures 

The acute need for high spectral resolution is highlighted by the 
detection of other species, which are producing signatures at al- 
titudes above the ozone level. 



5.3.1. Sodium detection 

The spectral region containing the Nal doublet is shown in 
Fig. [To] To repeat the discussion presented in Sect.[3l concerning 
the zero flux level correction, the result without any correction 
applied is also shown. Obviously, the atmospheric Nal doublet 
does not show up without that correction, and furthermore the 
broad wings of the two Nal solar lines are still present while our 
extraction process should have entirely washed out any solar line 
signature, as it is indeed almost perfectly the case in Figs.|6]or|9] 
This correction introduces an additional systematic error, 
which produces an altitude shift of about 1 .5 km per 0.0001 step 
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Figure 9. Observed binned variations (solid line) of the effective altitude h compared to model calculations (also binned). The 
Rayleigh-alone model calculation is shown (dashed line) along with the complete atmospheric model calculation (dotted line), 
which includes N2, O2 and O3 (see text). 



change in the zero level (see Fig.fTOll. This is small compared to 
the systematic error due to the instrumental order correction (on 
the order of ±2.5 km), which thus remains the major cause of 
uncertainty on our extracted altitude level. 

Note here a very interesting point. We know that both H2O 
and O2 produce strong absorption signatures in the Earth at- 
mosphere and, as shown in Fig. |2] they have relatively similar 
signatures near 5 900 A (in the Nal lines vicinity) for H2O or 
around 6 300 A for O2. We will see in the sections below that 
O2 is indeed easily detected in the grazing (transit like) observa- 
tions of the atmosphere, while we see here that none of the nu- 
merous H2O lines present in that spectral range (see Fig. |3]l are 
detected. The reason is very simple: H2O is hidden in the tran- 
sit like observations because it is at lower altitudes than ozone, 
while O2, a major constituent of the atmosphere, is still present 
in sufficient quantities at higher altitudes, and indeed shows up 
above the ozone layer. This is exactly what models do predict 
(see e.g.|Kaltenegger & Traub 2009 1 and thus reveals the quality 



and consistency of our analysis, from which only species thought 
to show up are indeed detected. 

The zero level correction induces additional very large errors 
on the estimated altitude of the Nal spikes: the Nal 5 890 A line 
peaks up from 35 to 100 km, while the 5 896 A one from 35 to 
65 km altitudes. 



Nal is known to be present in the Earth atmosphere within a 
layer at about 92 km, presenting an average thickness of about 
1 1 km jMoussaoui et al. 201 Oi l. If it was completely opaque in 
the two Nal lines, we should have found the same level in both 
lines at lower altitudes than the 92 km, depending on the ratio 
between the atmospheric line width (-0.03 A) and the instru- 
ment spectral resolution (-0.09 A). This should lead to reduced 
altitudes of about one third of 92 km above the 30 km ozone 
signature, i.e. ~50 km. 

This could correspond to our 0.0002 zero level correction, 
leading to both line levels at about 35 km altitude, a little low 
when compared to our evaluation. 

However, we can also directly evaluate the line's opacity 
in the 92 km altitude layer by using the observed average ver- 
tical sodium column density of the layer to be on the order 
of 4.10'"^ m"^ jMoussaoui et al. 2010 1. Using the average alti- 
tude, thickness, and "vertical" column density of the layer, we 
can translate these numbers into an average "horizontal" col- 
umn density at 92 km of altitude and found it to be on the or- 
der of 3.10" cm"^, leading to an optical thickness in the core 
of the line on the order of 3 in the strongest Nal line and 1.5 
in the weaker one. Both lines are thus opaque, but not strongly 
opaque. The known variability of the Nal layer could be the rea- 
son for this slightly less opacity, which renders it probable that 
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Figure 10. Nal doublet spectral region. The shifted h(A) variations (solid lines) are shown along with the level of the model ozone 
prediction (dotted line). Upper plot. The zero level correction as detailed in Sect.[3]is not applied. The two Nal narrow atmospheric 
signatures do not show up while through our complete data analysis method the broad wings of the sodium solar lines are obviously 
still present, contrary to expectation (the solar broad lines wings should not appear in any information uniquely related to the Earth 
atmosphere). Lower plot. The zero level correction by subtracting a 0.0003 value is now applied (see text). The Nal solar wings 
have disappeared, while the narrow Nal atmospheric signatures clearly show up. To show the impact of the zero level correction, 
the two evaluations corresponding to the 0.0002 and 0.0004 zero corrections are overplotted (thin dash-dotted lines). This gives a 
direct idea of the induced error bars on our evaluation. 



a marginally opaque layer was observed during the eclipse, then 
leading to two different altitudes for the observed doublet spikes 
corresponding to the also plausible values as the ones observed 
with our 0.0003 zero corrections, which are at 55 and 45 km 
respectively. 



From this brief quantitative discussion, we note that the 
0.0004 zero correction is certainly a limit because, in that case, 
the altitude of the strongest line is at more than 100 km, i.e. too 
high, even above the altitude of the layer itself. 

The Nal atmospheric detection is certainly compatible with 
what is known about the Earth layer, both qualitatively and quan- 
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titatively. This shows that our high-resolution approach could 
give us access to narrow line signatures, and thus gives us confi- 
dence in searching for other species. We detail this search in the 
following sections. 



Fig. |9]l, in which the model is higher by ~2 km. This band will 
certainly be scaled down by some amount in a lower tempera- 
ture model, simply because the atmospheric scale height will be 
reduced. 



5.3.2. The 6880 A oxygen molecular band 

Two molecular oxygen absorption bands have their signatures 
in the observed spectral range. We used for them the same zero 
level correction as for the Nal lines. 

The strongest one, the O2 forbidden '21^-^S„ transition band 

B (1-0) at 6 880 A, is shown in Fig. \n\ AU molecular band 
peaks are clearly present, undoubtedly signing the O2 detection. 
To evaluate the impact of the zero level adjustment, the three 
evaluations of h{A) made with the three values, 0.0003+0.0001, 
are also shown in Fig.[TT] This directly shows that in this spectral 
range, the zero level adjustment adds a systematic error of about 
~+2.5 km in the continuum and up to ~+10 km in the peaks, 
showing that their absolute estimates from these observations 
are quite imprecise. 

However, the average peak altitude and its relative height 
variations are certainly more meaningful. In particular, the av- 
erage height of the 0.0004 correction produces signatures closer 
to the model predictions, suggesting that the zero level correc- 
tion may be slightly wavelength-dependant or that its acceptable 
range is reduced to 0.00035+0.00005. One can also note that 
the observed peak heights when moving redward seem to drop 
more rapidly than in the model prediction, an obvious indica- 
tion of a cooler atmosphere than the "standard" one selected for 
the model calculation. This suggests that better adjusted Earth 
atmospheric model parameters in terms of vertical temperature- 
pressure profile are needed. This is however beyond the scope 
of this paper, in which we just intend to show that some species 
detection is feasible according to that observational approach in 
both qualitative and quantitative terms. The fine adjustment of 
model parameters will be done in future papers according to fu- 
ture observations to be made, all along a complete lunar eclipse 
to search for Earth atmospheric variations along the Earth limb 
sampled during the eclipse evolution. 

5.3.3. The 6 280 A oxygen molecular band 

The other oxygen molecular band in the observed spectral range 
is the y (2-0) at 6280 A. It is shown in Fig. [12] for which the 
same zero level corrections have again been used. 

The narrow spike positions of the band are clearly detected. 
This band is weaker however and emerges only by ~1 km above 
the continuum. A detailed discussion of this weaker signature is 
again beyond the scope of this paper, but this second detection 
of O2 demonstrates the quality and precision of our approach, 
which allows us to detect atmospheric signatures just beyond 
the photon noise. Indeed the photon noise induces in that region 
of high solar flux and high instrument sensitivity only a +0.5 km 
uncertainty over our altitude evaluation. 

Note that for this O2 molecular band we have checked the 
perfect match of the peak's spectral positions with the predicted 
ones while it seems that again their heights drop faster than 
the peak's altitude variations evaluated in the model calculation. 
This appears to confirm that the "standard" atmospheric model 
selected for that comparison is probably hotter than the observed 
atmosphere, particularly in the sampled ~30 km altitude range. 
If this were true it may also explain the slight mismatch found 
between the predicted and the observed O3 Chappuis band (see 
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Figure 12. y (2-0) 6 280 A O2 molecular band. Upper plot. 
Same background correction as in Fig. [10] (see text). The solid 
line represents the altitude information as extracted from the 
SOPHIE data set. Lower plot.The O2 band spectral structure 
is shown to demonstrate the clear coincidence in position of the 
repeated spikes signing the O2 detection. 



5.4. Previous lunar eclipse studies 

Other observations of the 2008 August 16 lunar eclipse were 
completed through inter-calibrated optical and near-infrared 
ground-based observations at the William Herschel and Nordic 
Optical Telescopes by Palle et al. (2009), to provide continuous 
wavelength coverage from 0.36 to 2.40 jum. Their approach was 
very similar to ours, except that their spectral resolution is some- 
what lower (~1 000 instead of ~75 000). They also focused on 
umbral, penumbral, and out-of-the-eclipse observations without 
detailing (as we did) the precise location within the penumbra 
from where the observations were made. This is important be- 
cause it leads toward only qualitative detections instead of quan- 
titative ones. 

Furthermore, the Palle et al. (2009) transmission spectrum 
is "calculated by computing the ratio of the umbra/penumbra 
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Figure 11. O2 forbidden 'ZJ-^Ej, transition band B (1-0) at 6 880 A. The same background correction as in Fig.fTOlhas been applied 
(see text). Upper plot. The solid line represents the altitude information extracted from the SOPHIE data, while the dotted line is 
the model calculation, in which O2 absorption is included from 6 856 A upwards. The narrow atmospheric signatures due to the O2 
absorption are clearly seen. The sharp peaks perfectly match the model with respect to spectral positions, while their heights could 
be quite different. Lower plot. The thick solid line again represents the evaluated altitude for the 0.0003 zero level correction, while 
the thin dash-dotted lines show the extreme possible variations of the estimated levels due to the 0.0002 (lower evaluations in the 
peaks and higher ones in the continuum) and 0.0004 (higher evaluations in the peaks and lower ones in the continuum) zero level 
corrections. 



regions taken at the same averaged air mass in order to mini- 
mize the local atmospheres telluric line variations". This means 
that the observed transmission spectrum they evaluated is not 
the transit-like spectrum of the Earth. Indeed the umbra spectra 
resulting uniquely from rays deflected by the lower and dense 
layers of the atmosphere are unobservable during a real transit 



of an extrasolar planet since the corresponding stellar light de- 
flected through the lower layers of the atmosphere moves away 
from the observer line of sight. What Palle et al. (2009) have 
reported are identifications of different atmospheric specie sig- 
natures only present in the lower parts of the atmosphere as seen 
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within the umbra, but not the real transit signatures, which are 
only observable within the penumbra, as we did. 




WAVELENGTH (A) 

Figure 13. The h(A) variation as extracted from the region of the 
Call H line. Because the S/N is quite poor in that region the 
photon noise in this part of the spectrum is more than 10 times 
higher than in the other spectral regions. Large signatures similar 
to the Nal ones should be visible however, but according to the 
low Call content of the Earth atmosphere, they are simply not 
observable, as it is the case here. 



In particular, Palle et al. (2009) have not detected the Nal sig- 
natures. They are produced at very high altitudes only (~92 km) 
and should be indeed unobservable in the lunar umbra. They 
mention the detection of Call signatures however, while it is also 
known that in the Earth atmosphere Call ions are only present 
at high altitudes dGranier et al. 1989] l. Call is furthermore about 
120 times fainter than Nal, thus the Call detection (and not the 
Nal one) is quite a surprise. Because we have in our hands much 
lower S/N SOPHIE observations at our disposal, we tried how- 
ever to look for possible Call signatures. As shown in Fig. [131 
nothing shows up. We conclude that the Palle et al. (2009) Call 
detection is probably an artifact. 

To enforce that argument we note that Palle et al. (2009) 
clearly detected H2O signatures (a low altitude species) 
at ~7200 A, while as shown from model predictions 



( Kaltenegger & Traub 2009j ), these could not show up in a tran- 
sit spectrum being hidden below the O3 and Rayleigh signatures 
(see Fig. |9]). This is why the different H2O signatures we have 
seen in our observation of the "vertical" Earth atmosphere (see 
Fig.|2]) do not show up in any of our "horizontal" transit-like sig- 
natures within the Earth atmosphere, because they are simply too 
low in altitude as predicted from the model calculations (see e.g. 
Fig. [To] where several H2O lines should have shown up around 
the Nal doublet, if present). 

The Palle et al. (2009) study in a larger spectral domain than 
ours reveals the great variety of detectable species, while ours re- 
veals how precisely the atmospheric content could be analyzed. 



5.5. Detection of O2 and O3 in an extrasoiar planet 
atmospliere during transits 

The purpose of this work is also to show and evaluate the fea- 
sibility of the detection of O2 and O3 in the atmosphere of an 
Earth-like extrasoiar planet during transits. 



We will not repeat the arguments of Lecavelier des Etangs 
& Ehrenreich (2005) here but only recall the extreme difficulty 
of these observations, which demand 10"^ - 10"** accuracy to 
be obtained over the observed stellar flux before, during, and 
after the planetary transit. This is a real challenge even merely 
in terms of photon noise, and asks for very large telescopes in 
the 10 to 30 m diameter range. This is only possible from the 
ground in a relatively short-term perspective, in particular in the 
frame of the ELTs (Extremely Large Telescopes) concepts that 
are planned to operate during the coming decade. 

The ELTs are however observing through the Earth atmo- 
sphere, and it is difficult to argue that such high precision accu- 
racies will be achieved. 

Our observations show that although the O3 Chappuis band 
seems to be easier to detect according to the broad spectral range 
over which it extends, it is still possible to think that uncontrolled 
systematics due to atmospheric fluctuations will perturb the de- 
tection. This is the obvious difficulty of broadband detections. 




VELOCITY SHIFT 



Figure 14. Correlation function extracted from the O2 forbidden 
'2+-^2j transition band B (1-0) at 6 880 A by using a mask 

centered on each of the spike positions over a ± 0.1 A width. 
The typical error bar of this function is indicated (arrow) and 
is on the order of the tracing thickness, according to the total 
spectral width sampled by the mask. 



On the other hand, our O2 observation via high-resolution 
spectroscopy signatures is very promising because it can be 
shown that the addition of all spectral lines over only the O2 
forbidden 'E+--'E^ ti-ansition band B (1-0) at 6 880 A is equiv- 
alent to a detection over a ~8 A band pass. This places the O2 
detection among the feasible ones from the ground with ELT 
telescopes at least in terms of photon fluxes. But the use of high 
spectral resolution and stable spectrographs, similar to the one 
presently used for radial velocity searches of extrasoiar planets, 
has proved to be possible in a quite difficult case, which is to 
observe the Earth's atmosphere as a transiting planet through the 
Earth atmosphere itself, i.e. all searched for and perturbing spec- 
tral signatures are exactly at the same wavelengths. For transit- 
ing extrasoiar planets the situation will be more favorable be- 
cause a simple shift of ~ 10 to 30 km/s or more will be enough 
to clearly separate them as shown through the correlation func- 
tion extracted with an O2 like mask (see Fig. [T4l) . This should 
be a relatively common observational situation. Note the sec- 
ondary peaks of the correlation function which are due to the rel- 
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atively regular separation between the successive O2 line peaks. 
Because these are however well separated, this should not pre- 
vent O2 detection over this molecular band. 

Our observation shows that with high spectral resolution, 
systematics are much easier to control because extremely nearby 
spectral regions behave similarly to each other in terms of any 
systematic effects. This may be one of the very few possible ap- 
proaches to reach this goal, which is detecting O2, one of the 
strongest signatures of life bearing atmospheres. 

Observing O2 and O3 species in extrasolar Earth like plane- 
tary atmospheres, both accessible from the ground in the same 
spectral range, shows that these searches will be attainable for 
ELT observatories equipped with high spectral resolution and 
high stabiUty spectrographs. 

6. Conclusion 

We have shown that Lunar eclipse spectral observations are able 
to reveal the atmospheric content of the Earth's atmosphere, in 
particular, broadband signatures of O3 and Rayleigh scattering 
as well as narrowband features of Nal and O2 observed at high 
spectral resolution. As these observations mimic future studies 
of transiting extrasolar planets, we are confident that quantitative 
information about extrasolar atmospheres will be within reach. 
Furthermore, both Rayleigh scattering and O3 are broadband and 
extend across the visible part of the spectrum, just where solar- 
like stars have their maximum flux, making it easier to obtain 
the high signal-to-noise values necessary for a positive exoplanet 
detection. 

More studies of Lunar eclipses should be completed in order 
to better quantify the present detections and more precisely show 
the feasibility of future extrasolar Earth-like planets studies. 
These studies, with extremely large ground-based telescopes, 
should allow at least the detections of ozone and O2. In any case, 
these observations will be extremely difficult because the re- 
quired accuracy is in the 10"^ to 10"^ range. This is why observ- 
ing in the visible range, where orders of magnitude more photons 
are available, could ultimately be one of the most promising ap- 
proaches. 
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